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Abstract 

It is shown that the axial and polar perturbations of the spherically symmetric black hole can 
be described in a gauge-invariant way. The reduced phase space describing gravitational waves 
outside of the horizon is described by the gauge-invariant quantities. Both degrees of freedom 
fulfill generalized scalar wave equation. For the axial degree of freedom the radial part of the 
equation corresponds to the Regge- Wheeler result |lj and for the polar one we get Zerilli result 
g], see also Q, E) for both. An important ingredient of the analysis is the concept of quasilocality 
which does duty for the separation of the angular variables in the usual approach. Moreover, there 
is no need to represent perturbations by normal modes (with time dependence exp(— iki)), we 
have fields in spacetime and the Cauchy problem for them is well defined outside of the horizon. 
The reduced symplectic structure explains the origin of the axial and polar invariants. It allows 
to introduce an energy and angular momentum for the gravitational waves which is invariant 
with respect to the gauge transformations. Both generators represent quadratic approximation of 
the ADM nonlinear formulae in terms of the perturbations of the Schwarzschild metric. We also 
discuss the boundary-initial value problem for the linearized Einstein equations on a Schwarzschild 
background outside of the horizon. 



1 Introduction 

In section 2 we introduce standard notions for linearized gravity and its 3+1 formulation. In the next 
section the analysis of constraints and Killing fields on an initial surface gives charges for the linear 
field on the Schwarzschild background. 

Section 4 contains the main technical results related to the so-called 2+1 decomposition of the 
initial data and to the gauge invariant description of the evolution. The invariants introduced in this 
section contain the full gauge-independent information about initial data. If we "insert" the initial 
value constraints into the canonical symplectic structure we can express the symplectic 2-form f2 in 
terms of the invariants. This is precisely shown in Appendix B. This way our invariants play a role 
of the reduced unconstrained initial data which is gauge independent. Moreover, the axial degree of 
freedom fulfills four-dimensional counterpart of the radial equation proposed by Regge and Wheeler 
@ > O > @ anc ^ t ne polar invariant is related to the Zerilli equation [|j , j| , Q . 

Section 5 is devoted to the stationary solutions, their behaviour on the horizon and precise inter- 
pretation of the "mono-dipole" solutions. In particular, the result of Vishveshwara || that stationary 
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axial perturbation can exist only for dipole perturbation can be easily recovered. However, if we as- 
sume more general conditions on a horizon we may have other solutions. The dipole part of the axial 
invariant corresponds to the spherical symmetry of the background metric and represents angular mo- 
mentum. On the other hand, there is no dipole invariants in the polar part because the background 
metric is not invariant with respect to the boosts and spatial translations^. Moreover, the lack of 
invariant in the dipole polar part corresponds to the fact that this part of the metric can be always 
"gauged away" . We would like to stress that mono-dipole perturbations of the Schwarzschild metric 
represent different phenomena than the higher-pole perturbations. It is clear from the approach that 

1 > 2 represents gravitational wave perturbation whereas I — and I = 1 correspond to the charges. 
For example, I = in polar part represents infinitesimal perturbation of the mass (we move in the 
space of solutions from one Schwarzschild solution to another), dipole part in axial part corresponds 
to the infinitesimal angular momentum (we move from Schwarzschild to Kerr solution) . Formally also 
monopole part in axial degree of freedom can be analyzed and represents Taub-NUT charge but this 
move is excluded by topology. 

Reduction of the symplectic form presented in Appendix B allows to introduce invariants from 
the symplectic point of view in section 6. In the next section we define (in a gauge-invariant way) 
hamiltonian systemF], energy and angular momentum generators and boundary-initial value problem 
for the linearized field on the Schwarzschild background. Moreover, we are able to show that the 
obtained hamiltonian is a quadratic approximation of the ADM mass defined at spatial infinity for 
the full nonlinear Einstein theory and this particular result will be described elsewhere. 

The energy and angular momentum generators are well defined for regular radiation data x, 2L Yj 
Y which is finite on the horizon and vanishing at spatial infinity according to the S.A.F. Christodolou- 
Kleinerman condition IQ. 

In Appendix A we show how to reconstruct the full four-metric from the invariants assuming 
the gauge conditions used in and Q . This construction explains the relations between our invariants 
and the special form of the metric used in jjj and jj). In particular, we examine precisely "mono- 
dipole" part of the metric (1 = and I = 1) which seems to be not fully analyzed in literature. 

2 Linearized gravity 

In this section we remind some standard notions related to the Einstein equations and the initial value 
problem. 

Linearized Einstein theory (see e.g. [0 or ||) can be formulated as follows. Einstein equation 
2G^(g) = 16^V (2.1) 
after linearization gives 

W ;Q + h ^' a - - (v aP h af} ).^ - ri ia ,[h a p' af> - h a a '%] = 167rT M „ (2.2) 

where pseudoriemannian metric g^„ = rj^ + h^v, r\^ v is the background metric and ";" denotes four- 
dimensional covariant derivative with respect to the metric rj^ v - Moreover, we assumed that 7^,„ is a 
vacuum solution of Einstein equation (G M „(?7) = 0). 

The (3+l)-decomposition of (|]^) gives 6 dynamical equations for the space-space components hki 
of the metric (latin indices run from 1 to 3) and 4 equations which do not contain time derivatives of 
hki- It is possible to introduce this decomposition straightforward. However, a natural way for the 



formulation of the (3+l)-splitting for the equation (2.2) is to linearize ADM formulation of the initial 
value problem for the full nonlinear Einstein equation (^j]). We shall introduce canonical variables 
for the linearized case. They appear in a natural way if we start from the ADM formulation of the 
initial value problem for Einstein equations . 

1 Pin ole invariants exist in flat Minkowski space and they represent linear momentum and center of mass respectively 

(sec Fa). 

2 After separation of the angular variables several "hamiltonian" results proposed by the author can be easily trans- 
lated into Moncrief's approach presented in [H. 



2 



Let (gki,P kl ) be the Cauchy data for Einstein equations on a three-dimensional space-like surface 
Ti. This means that gki is a Riemannian metric on £ and P kl is a symmetric tensor density, which we 
identify with the ADM momentum |Q, i.e. 



P M = Vdetg mn (g kl TrK - K kl ) 

where K k \ is the second fundamental form (external curvature) of the imbedding of £ into a spacetime 
M. 

The 12 functions (gki,P kl ) must fulfill 4 Gauss-Codazzi constraints 



Pi l \l = 87rv/det g mn Ti^ (2.3) 

(detg m „)R- P kl P M + \(P kl 9ki) 2 = l6ir(det Smn )T MV nf i n v (2.4) 

where T^„ is an energy momentum tensor of the matter, by R we denote the (three-dimensional) scalar 
curvature of gki , is a future timelike four-vector normal to the hypersurface £ and the calculations 
have been made with respect to the three-metric gki ("|" denotes the covariant derivative, indices are 
raised and lowered etc.). 

The Einstein equations and the definition of the metric connection imply the first order (in time) 
differential equations for gki and P kl (see [[z] or M p. 525) and contain the lapse function N and the 
shift vector N k as parameters 

2N ( 1 



IH = -j= ( Pkl - 2&kiP ) + N k{l + N l{k (2.5) 
where g := detg m „ and P := P kl g ki 



P kL = -N^/gR kl + ^g [N lkl - g fc 'iV |m 



~AVgg fe 'R- ^ (p km P m l - ipP W ) + {P kl N m ) lm + 

A ' ' f P kl P k i - \p 2 \ - N k Wl P ml - N l ]m P mh + ^NJgT mn g km g ln (2.6) 



2VI V 2 

Let us consider an initial value problem for the linearized Einstein equations on Schwarzschild 
background rj^: 

2m \ , o / 2m x 1 



r 



rj^dx^dx" = - 1 dr + 1 dr A + r 2 d9 2 + r 2 sin 2 6d<p 2 (2.7) 



together with the radial coordinates: x 3 — r, x 1 = 9, x 2 — ip. Moreover, t = x° denotes the time 
coordinate. We consider only the part of the Schwarzschild spacetime outside of the horizon, r > 2m. 

We use the following convention for indices: greek indices /i, v, . . . run from to 3; k,l, . . . are 
spatial coordinates and run from 1 to 3; A, B, . . . are spherical angles (9, (p) on a two-dimensional 
sphere S(r) := {r = x — const} and run from 1 to 2. Moreover, let T)ab denote a two-dimensional 
metric on S(r). 

Let v := 1 . There are the following non- vanishing Christoffel symbols for the metric (|2.7|): 



.3 m . r 3 v A 1 A 3 tyiv q m A 

33 — n j t AB — VAB , 1 3B = -0 S , t 00 = T" j 1 30 = n , * BC 



where 5 A b is the Kronecker's symbol and T A bc are the same as for a standard unit sphere 5(1) (in 
usual spherical coordinates T 6 \<t> = — sin cos and P^^g = cot 9). 
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The curvature of the background metric we denote by C p V \ K and the following components of the 
Riemann tensor are non- vanishing (up to the symmetries of the indices) 



CO r~iZ 
AOB — O A3B 



in 



c 



AB 



CD 



2m , j 



s;A cB 
70 ' 



5 A D 5 B c) 



We can introduce the following submanifolds of the Schwarzschild spacctime M: 
x° = s, x 3 > 2m} - 



{x G M 



\J S s (r) where S s (r) := {x G S s 

r£[2m,oc[ 



r} (2.8 



and S s is a partial Cauchy surface outside of the horizon. 
The surface £ carries the induced Riemannian metric rjki ■ 

1 



i]kidy k dy l 



(2.9) 



Usually it is convenient to change the coordinate r to which is a solution of an ordinary differential 
equation (see Q): 



r + 2m In f 1 ) 

V 2m / 



(2.10) 



dr* 1 
-T- = - ; r, 

dr u 

and moves horizon to — oo. 

The ADM momentum P kl for the metric (2.9) on each slice S s vanishes (P fc ' = 0). The shift vector 
is also trivial (N k = 0) and the lapse N = ^/v is vanishing on the horizon. Moreover, the Ricci tensor 
for the three- metric rju has the following components: 



Rfcz — 



1 



- r N\ 



\kl 



and the scalar curvature R vanishes. 

Let us define the linearized variations (hki, P kl ) of the full nonlinear Cauchy data (gki,P kl ) around 
background data (rjki, 0) 

hki ■= gu - Vki , P kl ■= P kl 



(2.11) 



We should now rewrite equations ([2.3])-(2.6) in a linearized form in terms of (hki,P kl ). Let us 
denote P := rjkiP kl and h := rj kl hki- The vector constraint (|2.3|) can be linearized as follows 



Pi l \i 



P 



' I' 



(2.12) 



Let us stress that the symbol "|" has different meanings on the left-hand side and on the right-hand 
side of the above formula. It denotes the covariant derivative with respect to the full nonlinear metric 
ghi when applied to the ADM momentum P kl , but on the right-hand side it means the covariant 
derivative with respect to the background metric rjki- The scalar constraint (2.4) after linearization 
takes the form 



1 



VgR-^ (P fe 'P fe * 



1 



Vv 



h kl R k i 



where y/rj := ^/detrjki. 

The linearized constraints for the vacuum (T^ 



0) have the form 







h kl R k i = 



The linearization of (|2.5D leads to the equation 



2N ( 1 



l 0k\l 



l 0l\k 



(2.13) 

(2.14) 
(2.15) 

(2.16) 



4 



where N :— 



= ^/v, N k = rjok = are the lapse and shift for the background. Let us denote 



the linearized lapse by n := ^y/vh®. The linearization of (2.6) takes the form 



-nRki - NSRki + n\ u - N m 5T r ' 



kl 



'iVi 



\mn 



(2.17) 



where 



ST r > 



kl 



k\i 



+ h r ' 



i\k 



is the linearized Christoffel symbol and similarly 



SR 



kl 



2 {h m k\lrn 



1 1 km 



*'kl\m 



- h 



\kl) 



is the linearized Ricci tensor. 

It is well known (see for example [|l9|) that the linearized Einstein equations (2J2) are invariant 
with respect to the "gauge" transformation: 



(2.18) 



where ^ M is a covector field. The (3+l)-decomposition of the gauge acts on the Cauchy data in the 
following way 

2vA- 1 P k i -» 2vA- 1 Pki + (vtf k )\i + {v§) lk -Vki« olm ) m (2.19) 

hki -> fc*l + 6|fc + &|i ( 2 - 2 °) 

where A := ^/v^/dct r\ki (= r 2 sin0). It can be easily checked that the scalar cons train t ( 2.15|) an d 
the vector constraint ( 2.14 ) are invariant with respect to the gauge transformations ( 2. IS ) and ( 2.20 ). 

The Cauchy data (hki, P kl ) and (hki, P ) on E are equivalent to each other if they can be related 
by the gauge transformation £ M . The evolution of canonical variables P kl and hki given by equations 
( 2.16 ), ( 2.17 ) is not unique unless the lapse function n and the shift vector hok are specified. 

We will show in the sequel that it is possible to define a reduced dynamics in terms of invariants, 
which is no longer sensitive on gauge conditions. The construction is analogous to the analysis given 
in Hi. 



3 "Charges" on the Schwarzschild background 



The vector constraint (2.14) allows to introduce "charges" related to the symmetries of the background 
metric. There are three generators of the rotation group, which are simultaneously Killing vectors on 
the initial surface S. Let us denote this Killing field by Z k . It is a solution of the Killing equation 



Zk\i + z, 



l\k 







(3.1) 



Let V C £ be a compact region in E. For example V := S s (r) and dV — S s (ro) U S s (r\). 

r£[r ,ri] 



From ( |2.14[ ) and ( J37l| ) we get 

0= f P k \iZ k = [ (P kl Z k \^ [ P 3k Z k 



(3.2) 



The equation (3.2) expresses the "Gauss" law for the angular momentum charge "measured" by the 
flux integral. It is easy to relate this charge to the dipole part of invariant y, which will appear in the 
sequel - (4.19) in the next section. For example, when Z = d/d<p we have 



16tts 2 



:= 167rj^ = -2 [ P\ = -2 / P\(r 2 e AB cos e\ B 

JdV JdV 
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2/ r 2 P 3 AllB e AB cos9 = / Aycosfl 
Jov Jav 



(3.3) 



The time translation defines a mass charge from the scalar constraint (2.15) as follows 



(fc W |,-fcl fc ) 



h kl N\ 



\kl 



NVV (h lk \k ~ fc' 1 ) + ^ (N\ l h - N ]k h H ) 



av 



\ I h 3k {k - h 13 + ^hN^ 3 - ^h k3 N v ,. 



N 



N 



(3.4) 



and it can be related to the monopole part of an invariant x introduced in the next section by (4.21). 
16ttp° = 



av 



A (Nh 3k \ k ~Nh} 3 + N k h k3 - N 3 h"j 

3m, \ 
vrH, 3 -(l )H) = 



- 2h 33 
av t \ 



A 



B _1 x 



av 



(3.5) 



4 Equations of motion for the invariants 



The radial foliation of the Cauchy surface £ related to the spherical symmetry allows to perform 
(2+l)-decomposition of the initial data. In this section we introduce reduced gauge invariant data on 
£ for the gravitational field, similar to the invariants introduced in [fL6f . For this purpose we use a 
spherical foliation of £ (see formulae ( 2^8 ) and (2.S)). 

o 

Let A denotes the two-dimensional Laplace-Beltrami operator on a unit sphere S(l). Moreover, 
H := r) AB h A B, Xab ■= ^ab-^abH, S := rj AB P A B, Sab ■= Pab-^VabS according to the notation 
used in fig] . 

The spatial gauge ( [2.2C| ) splits in the following way 



1*33 -> h 3 3 + — ^{VV&),3 

V v 

2 

h-3A — > h,3A + &,A + £a,3 £,A 



h-AB -> h AB + £a\\b + €b\\A + -Vab£, 



(4.1) 
(4.2) 
(4.3) 



where by "| |" we denote the covariant derivative with respect to the two-metric t\ab on S(r). Similarly, 
the temporal gauge (2. IE) can be splitted as follows 

(4.4) 
(4.5) 

(4.6) 
(4.7) 



A 


-1 P 3 3 


-> A - 


1 P 3 3 


-e 0|l V 


_ £ ,0,3 

r 


A 


~ X P3A ~ 


+ A" 


l P3A 


+ (# - 












r vr l J M 


A" 


^Sab - 


■* A~ 


1 Sab 




~ \VAB^ C C 




Kr l S - 


■* A" 


1 S- 


~(Ne> 3 


),3--C 0,3 -C 0||c 

r 



It is also quite easy to rewrite the (2+l)-decomposition of (2.16) 



, 33 =A- 1 (F 3 3-5) + — (^ 03 ),3 



(4.8) 
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-1 * 
h 3A = 2vA P 3 A + ho3\\A + h A,3- ~h A 



h A B = 2vA Sab-VabA P ' + <W||b + Kb\\a + 2w r} AB ho3 
The dynamical equations ( [2.17 ) take the following (2+l)-form: 



(4.9) 
(4.10) 



V V 



(H, 3 - 2h, A u - 2r- 1 /i 3 J ), 3 + 2r- 1 (ff, 3 - 2^p - 2^^/) 



(4.11) 



2A- 1 Ao = 



^(^°o),3-r- 1 /i°o 



lie 



-^^ ||c + -(iJ,3-2/,3 A ||A-2r-^ J 3 ) 



+ ^3C||yl"' 4 + — ^3C - X A C||A, 3 



(4.12) 



2A ^ Sab = h° oUB - -^ab^o 11 ^ + h 3 3 \\AB - \vA B h 3 3^ C c+ 

-{h 3 A\\B + h 3 B \\A -VA B h 3C \\ C ),3 + {vX°B,3VCa),3 + 

IIC C C CD ^ 

+ XAB ||C* — X a\\bc - X B\\AC + VABX \\CD + ^XAB 

2A" 1 5 = -X^, 3 ),3 - ^,3 - + (^3 + ^) i|A A + ^(^3 + H) ^h 3 

2m 



(4.13) 



-If, 3^ 



+ -5-^3,3 + 3 - 2h 3 A \\A - 2r- i / l3 d ), 3 + 4-(if, 3 - 2/i 3 ||A - 2r- 1 ft. 



2X 



CD 



1 1 CO 



(4.14) 



The vector constraint (2.14) splits in a similar way 
-^(V^P 3 3),3 + P3 A l \A-r- 1 S = 



P 3 a, 3 + Sa b \\b + 2^W A — ^ 



(4.15) 



(4.16) 



And finally the (2+l)-decomposition of the scalar constraint (2.15) can be written in the form 



\kl 



<3 A 
1 3||A 



2r- 2 h 3 , 



6m 



h 3 3 



MrH,3~2rh3A nA ~2h 3 3 ) 
1 



,3 + 



H^ c c + r- 2 H - X AB \\AB = 



(4.17) 



Let us notice that we can split the dynamics into two separate parts which we call axial and 
polar respectively. The axial part of initial data consists of two momentum components P 3A ^ b eab, 
S c a\\cb£ AB and two metric components h 3 A\\B^ AB , X° a\\cb£ AB ■ The only gauge freedom is con- 
tained in £,a\\b £AB which acts on the metric components. The gauge invariants P 3A ^ B eab and 
S c A\\CB eAB are related by the curl part of the vector constraint (4.16) as follows: 



(r 2 P 3A ^eAB),3 + r 2 S C AHCBS AB = 
It is easy to verify that the following pair of gauge invariants: 
y := 2K- 1 r 2 P 3A W B e A B 



(4.18) 
(4.19) 
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(4.20) 



contains the whole information about axial part of initial data up to the gauge freedom (see also 
Appendix A). We will show in the sequel that this is a canonical pair with respect to the symplectic 
structure of linearized Einstein equations. 



Let us define the polar invariants as follows 
x := r 2 X AB UB - ~( A + 2)H + B [2h 33 - 



2r 2 S AB UB + B 



2rP 



3A 



AP 3 3 



2rh 



3C 



lie 



rvH, ; 



(4.21) 
(4.22) 



where 

° . ( ° 6m\ 1 

B:=(A+2) A+2-- 

is a quasilocal operator - it is local with respect to the coordinates t, r but non-local on each sphere 
S(r). The proof that x and X contain the whole information about polar part of initial data P 3A \\Ai 
S AB \\ab, -P 3 3, S, h 3C \\ c , X AB \\ab, h 3 3 , H up to the gauge transformation £°,£ 3 ,£ A ||a we give in 
Appendix A. Moreover, in Appendix B we show that (x, X, y, Y) is the reduced canonical initial data 
on S. 

We can check the reduced field equations for the axial invariants 



1 Y 



A 



r(r 2 y), 



, 3 + ^(A+2)y 



(4.23) 
(4.24) 



More precisely, the curl part of (4.12) gives (4.23) and (4.24) follows directly from (4.9), (4.10) and 
may be checked by inspection. 

It can be easily verified that the axial invariant y fulfills the generalized scalar wave equation (as 
a consequence of the above dynamical equations) 



□ + ^)y = o 



(4.25) 



where D denotes the usual wave operator with respect to the background metric r]^ u . 

There exists a simple relation between the equation (4.25) and the so-called Regge- Wheeler equation 
@i 01; |1 ■ Let us rewrite equation (4.25) in the following form: 



-y + -Nry), 3 ]. 3 = ^ ( } y 



(4.26) 



where ' is a "spherical operator" defined as follows: 



v , 
■-t A 



6 m 
r 



If we insert the invariant y in a special form y = exp(iat)Yi(9, <f>)Z^ '{r)/r into equation (4.26), we 
obtain Regge- Wheeler equation 



dr 2 



Here Yj is a spherical harmonics such that A Y\ = —1(1 + l)Yj and coordinate r* is defined by (2.10). 
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The polar invariants fulfill the following equations: 



-X 
A X 



A 



{( 



w 2 x, 3 ) 



A +v(l - 2B) + 1 6x| 



(4.27) 



(4.28) 



To ob tain ( [4.27[ ) we need e quations (p~g| ), (p~9|), ( |4. ldj ) and v ector constraints ( |4.15| ), ( [4.f6[ ). Similarly, 
(4.28) is a consequence of ( |hll| - fhl~3 ) and scalar constraint (4.17) (see also Appendix A). 

There exists also a generalized scalar wave equation for the polar invariant but it is no longer local, 
it is only quasilocal: 



i—i 8m. ° 
□ + -5- A 



A 



r 



6m 
r 



x = 



The similar equation to the ( 4.26 ) can be presented in the analogous form 



X+-K rX ),3], 3 



but now the operator y(+) is defined "quasilocally" 

6m \ 36m 2 



(4.29) 



(4.30) 



v 



( A + 2) A 



A +2 



2m 
r 



A 



6m 
r 



exp(icrt)Yi(8, <f))Z^ + \r) /r into equation 



+ a 2 Z (+) =^ (+) ^ (+) 



If we insert the invariant x in an analogous special form x 
( |4.30D we obtain Zerilli equation ||, |], g 

dr 2 

This way we have shown that both equations (Regge- Wheeler and Zerilli) posses gauge invariant 
formulations and their solutions contain the entire gauge-independent information about the linearized 
gravitational field on the Schwarzschild background (see also @). 

Let us notice that x and y are scalars on each sphere St{r) with respect to the coordinates x A . 
For the scalar / on a sphere we can define a "monopole" part mon(/) and a "dipole" part dip(/) as a 
corresponding component with respect to the spherical harmonics on S 2 . Similarly, the "dipole" part 
of a vector v A corresponds to the dipole harmonics for the scalars v \\a an d £ Va\\b- Let us denote 
by / the "mono-dipole-free" part of /. According to this decomposition we have 



y = dip(y) +y 

The dipole part of x and monopole part of y are vanishing ^. T he rest of the mono-dipole part of each 
scalar can be solved explicitly from the equations ( 4.23 )-( 4.28 ) and the solution has the form 



4m 

r — 3m 



12s 



y = 



From (4.23), (4.27) and the observation that mon(X) = dip(Y) = we obtain 







This is included in the definitions (4.21) and (4.19). More precisely, y is a divergence (see also (B.l)) and 



r 2 X AB \\AB + { A 



+ 2) ...]. Moreover, dip(x AE ' \\Ab) = because double-divergence of any traceless tensor is 



"mono-dipole-free" (see also (B.2) in Appendix B) 
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Moreover, Am = 0, ( A + 2)s = 0, which simply means that m is a monopole and s is a dipole, and 
they are constant with respect to the coordinates t, r. They correspond to the charges introduced in 



section 3. More precisely, m = p . Moreover, the angular momentum charge (3.2) can be obtained 



from the relation between spatial constant three-vector in cartesian coordinates and dipole harmonics 



s l zi 

s = 

r 



where zi are cartesian coordinates, r — \J S kl ZkZi and s l is a corresponding three-vector representing 
angular momentum (see Jh|). 

5 Stationary solutions 

For the axial degree of freedom y we can rewrite equation ( 4.26| ) using a new coordinate z := 221. 



4rn 2 y = (1 - z) 2 z 4 ^ - (1 - z)z i ^- + (1 - z)z 2 ( A + 4z)y 
oz z az 

This way horizon corresponds to z = 1 and spatial infinity to z = 0. Let us consider stationary 
solutions of the above equation which are regular at the spatial infinity (corresponding to z = 0). 

(l-z)z 2 ^-z 2 ^ + ( A+4z)y = (5.1) 
oz z az 

Moreover, if we separate the angular variables and include standard asymptotic behaviour at z = 0: 

y = z l+1 Yi{x A )u{z) 
then the equation for the function u is relatively simple 

(1 - z)zu" + [21 + 1 - (21 + 3)z]u' -(I- 1)(/ + 3)u = 
and the solution regular at z = is given by the hypergeometric function 

u = F(l-l,l + 3,2l + 2;z) 
In particular for I — 1 the function u is constant, it represents the angular momentum charge solution 



(3.3), and corresponding y is finite on the horizon. On the other hand, for I > 2 we obtain logarithmic 



divergence of the hypergeometric function F at z = 1. More precisely, 

F(l - 1, 1 + 3, 21 + 2; z) = z^ 2 '" 1 [P{z) + P(z) ln(l - z)\ 

where P and P are polynomials. The solution is not regular at z = 1 and it confirms the result of 
Vishveshwara that the only nontrivial stationary perturbation is given by the axial perturbation 
with 1 = 1. The interpretation of this solution is given by ( |3.3| ). One can check by direct computation 
for the Kerr metric [[To) , |nj : 

p 1 := r 2 + a 2 cos 2 8 A := r 2 - 2mr + a 2 

2mr 2mra sin 2 9 p 2 2 

.9oo = -1 + — jj- ffo0 = 9rr = gee = P 

■ 2 /> ( i 2 2mra 2 sin 2 6* 
= sm" I r + a H 

that the infinitesimal angular momentum gives the same result. More precisely, the linear part of the 
Kerr metric g^ v with respect to the parameter o: 

g^dx^dx" = rj^dx^dx" - 4— sin 2 9dtdcf> + 0(a 2 ) 

r 
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gives only the one non- vanishing component of the linearized Kerr metric on the Schwarzschild back- 
ground: 



h-Ocf, — 



2ma 



■ sm 



From the definition (4.19) and the equation ( |4.9| ) we can calculate the invariant y = 12ma cos 6*/r. 
It is easy to compare the result with (3.3) and it gives angular momentum charge s z = ma. 

Can we consider monopole solution in axial invariant? The definition (4.19) in terms of the initial 
data does not allow nontrivial monopole part of the axial invariant. However, we can consider such 
situation if we admit singular metric. Formally this case corresponds to the infinitesimal Taub-NUT 
solution. The Taub-Nut metric [O, fLl: 



g^dx^dx" = (r 2 + P)(d8 2 + sin 2 9d(j) 2 ) + {PMr 2 - v(dt + 21 cosfld^) 2 = 

v ~- 2 /77 r — I 
= ij^dx^dx" - Avl cos Qdtdcj) + 0(l 2 ) ; v := - 



Al 3m 

gives the linearized metric hos = —2lvcos0 and the monopole axial invariant y = (1 ). We 

r r 

should stress that this is only formal calculation because the tensor = — 2lv cos 9 is not well defined 
along the z— axis and is excluded as a global solution. The monopole charge in y plays a role of the 
topological obstruction for the existence of the global metric. It is similar to the magnetic monopole 
in electrodynamics (see also |i~9|). 

For polar degree of freedom x from equation ( 4.30| ) we get 



4m 2 x = (1 - zfz A -^ - (1 - z)z i ^ + (1 
oz z oz 



1 2 , 

3 + 3<A-l)(A 



■ 2)( A + 2 — 3z) 



( A + 2)x 



and the corresponding stationary equation has the following form 



(1 _ z )z 2 — - z 2 — 
1 ' dz* dz 



-(A-l)(A+2)(A+2-3z)- 2 



( A + 2)x = 



(5.2) 



For I = the solution x = 



precisely, if we put in the metric 
we get 



is related to ( |3.5| ) and corresponds to the mass charge. More 
m + Sm instead of parameter m and take the linear part in Sm 



2Sr 



2Sr 



Vn»( m + Sm)dx^dx u = ■q iiu {m)dx il dx u H dr H — dr 2 + 0(Sm 2 ) 



and the invariant x = 2Bh 



ASr 



3m 



r 

-l 



1 ) . If we compare with mass charge (3.5) we obtain 



p° = Sm (see also pj). 

For I > 2 we have the following transformation law: 

o o f)v r o o o 

A( A + 2)x = -6z 2 (l - z)^- + A(A+2) + 6z(l - z - 18z 2 (l - z)( A + 2 - 3zY 

oz L 



which moves the solution of (5.1) into solutions of (5.2). In particular, it is clear that for / > 2 the 
stationary solutions of the equation (|]^) have the same logarithmic divergence on the horizon as the 
solutions of the equation (5.1). The explicit stationary solutions (in a specific gauge) were also given 
by Zerilli in |g{ and here we present a gauge invariant confirmation of his result. 

Remark Although I = 1 is excluded in the definition of x, we can consider another variable 



(A 



(A+2)-VY B p B 



1 

2 H 



° 6m 
A +2 

r 



12^ 



2rh 



3A 



rvH ■ 
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which is no longer gauge- invariant in its dipole part (see also (A..8) in Appendix A). More precisely, 
it transforms with respect to the gauge transformation as follows 

dip (^H + ^-vQ^j = dip(-( A + 2)- 1 x) — ► dip(-( A + 2)- 1 x) + dip(£ A ,| A ) 

where Q is defined in Appendix A by ( |A.2| ). Formally, the dipole solution ( A +2) -1 x = ln(l— z)Yi{x A ) 
fulfills the same equation fl5.2| ), and it corresponds to the Regge- Wheeler infinitesimal translation as 
a gauge transformation^ (eq. 32 in 

„ sin6> / 2m\ ^ A cos6 , /„ 2m 

e = cos6; ^ = In 1 ; C A \\A = In l 

Ira \ r J m \ r 

This way the stationary polar dipole solution corresponds to the infinitesimal translation gauge and 
it is also logarithmically divergent on the horizon. 



6 The symplectic form and its reduction 

In this section we show the relation between the symplectic structure and the invariants introduced in 
the section 4. Let (P kl , hu) be the Cauchy data on a hypersurface S. Let us consider the symplectic 
structure 

fl := I 5P kl A5h k i 

It is invariant with respect to the spatial gauge transformation which is fixed on the boundary 
(<^las = °) : 



5P kl A Sh kl / SP kl A 6h kl +2 SP 3 t A 5? 

is JOT, 

Moreover, it is invariant with respect to the temporal gauge if we assume that £° and its normal 
derivative^] are fixed on the boundary: 

SP kl A 8h kl f SP kl A 5h kl + 

is 

+ [ VV U( N f)\k A 8h 3k - S(Nff A 8h + N5£° A S(h^ 3 - h 3 \i) 

Roughly speaking, the symplectic structure is invariant with respect to the gauge modulo boundary 
terms. The quadratic form J s SP kl A Sh k i can be decomposed into monopole part, dipole part and 
the remainder in a natural way. 

From the considerations given in the Appendix B (formulae ( p3.5| ) and ( B . 1 3| ) ) we can easily see 
that the "radiation" part 

Q = I SP_ kl A 5h kl ~ f 5XA A~ 1 ( A + 2)~ 1 5x + 5Y_A A _1 ( A + 2)~ 1 <5y 

JS JS 

where symbol "~" denotes equality modulo boundary term. Moreover, the "mono-dipole" part has 
the form (see (B.E) and (B.4) respectively) 

mon( / 5P kl A 8h kl ) = [ ^6P 33 A <56" 1 mon(x) + \ I rSP 3 3 A Smon(H) (6.1) 
Js is 2 2 J gs 



4 If we assume that the translation corresponds to the £ 3 = cos 8 then the component £ 9 is uniquely defined jas the 
polar gauge transformation preserving the gauge condition hs A \\A = which has been used by Regge- Wheeler and 



Vishveshwara 

5 You may "correct" the symplectic form by a boundary term in such a way that the result is gauge invariant for ■ 
vanishing on the <9£ without an extra assumption about derivatives (sec |20|). 
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dip( f SP kl A Sh kl ) ~ - / A<5dip(y) A A~ 1 6(h 3A \\ B e AB ) (6.2) 

The monopole part of the polar invariant mon(x) and the dipole part of the axial one dip(y) 
represent mass and angular momentum respectively, and they are supposed to be fixed. These quan- 
tities are analogous to the electric charge in electrodynamics. If we assume that there is no matter 
inside volume £ then both of them are fixed by the constraints, provided that they are controlled 
at r = 2m as a boundary condition on the horizon. Let us assume that <5rnon(x)| iS / r _ 2m ^ = 0, 
<^dip(y)| S / r _ 2m ) = then the mono-dipole part (6.1) and ( |6.2| ) vanish and the symplectic structure 



reduces to the "mono-dipole-free" invariants 

f 5P kl /\5h kl ~ [ <5X A A _1 ( A + 2)-~ L 5x + 5Y A A _1 (A+2)-^y (6.3) 

This way we obtain X, x, Y, y as the quasi-local canonical variables describing reduced unconstrained 
initial data on S. 

Remark One can ask the question when ( |6.3| ) is a strict equality not only modulo boundary 
term. The symplectic 2-form ft reduces to the mono-dipole-free invariants if we assume the following 
boundary conditions which fix the gauge freedom on the boundary: 

6h AB \ 9 z = 8{2h\ + 2rhz A u - rH l3 )\ = 



8{P iA ^e AB ) =0 5{Pl + 2rk 1 p 3A \\ A , 

3E " 



= (6.4) 

as 

In Appendix A we analyze the possibility when H_, Q, II and X AC \\cb^a B are precisely the gauge 
conditions and then X \\AB corresponds to the control of x. Roughly speaking, the "control mode" 



given by (3.4) contains four boundary conditions related to the gauge freedom plus two boundary 
conditions for the unconstrained degrees of freedom which we propose to call Dirichlet boundary 
data. 

If we introduce quasilocal coordinates: 

q 1 := A- 1/2 ( A + 2)" 1 / 2 x; p x := A~ 1/2 ( A + 2)~ l / 2 X 
q 2 := A- 1/2 (A+2)- 1 /2y. p2: = A -V2 (a+2 )-V2y 
we can rewrite the reduced symplectic structure ( |6.3| ) in the canonical form 

2 

ft = / SXA A _1 ( A +2)- 1 Sx + SY A A _1 ( A +2)- 1 <5y = V / §p n A Sq n 



Unconstrained initial data for the full nonlinear theory (which is similar to the considerations in this 
article) has been proposed in |l8). Moreover, the concept of quasilocality appeared in JTtJ and has 
been developed in pof . The boundary data possesses its counterpart in the full nonlinear theory (see 
proposition in [^o| ) and it will be discussed elsewhere. 

7 Energy and angular momentum of the gravitational waves 



The reduction of the symplectic form (6.3) allows to formulate the hamiltonian relation in terms of 
the reduced canonical variables 



/ x A _1 ( A + 2)- 1 JX - X A _1 ( A + 2)-Mx+ 
+ [ yA~ x ( A +2)-!^- Y A _1 ( A+2)- 1 6y = 16ir5H+ 
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/ -«(rx), 3 A-\ A +2)- 1 <5x+- U (ry), 3 A _1 ( A + 2)"% (7.1) 



where 



1 f v i ° i i 

16ttW := - J — X A ( A + 2)- x X + -Y A _1 ( A + 2)- J Y+ 

2 

«(rx), 3 A _1 ( A +2)- 1 (rxX 3 + x— V {+) A~ x ( A + 2) _1 x 



1 r A 



2 7 s r 2 
- [ — 



v 



v(ry), 3 A _1 ( A + 2)- 1 (ry), 3 + y-V^ A _1 ( A + 2) 

v 



-v 



(7.2) 



(see also eq. 4.19 and 5.34 in [£§). 

Similarly for angular momentum we propose the following expression 



16wJ(Z) = [ XA- 1 (A+2)- 1 Z A d A x + YA- 1 (A+2)- 1 Z A d A y (7.3) 



where Z is a Killing field ( |3.l| ). In particular for Z = d/d(f> the z— component of the angular momentum 
takes the form: 

16ttJ z = /xA- 1 (A+2)- 1 x. + YA- 1 (A+2)- 1 y,0 (7.4) 

The conservation laws for the energy and angular momentum 
d H = 
d J(Z) = 

are fulfilled if we assume appropriate boundary conditions on the horizon and at the spatial infinity. 
The natural choice of those conditions is to assume that x and y are vanishing on the boundary dT,. 

After separation of the angular variables in ( |7.2| ) we obtain a hamiltonian which has been used for 
the energy method by Wald H (see also in Q eq. 386) and it confirms stability of the Schwarzschild 
metric (see also ||). Here we have shown how to combine an energy of different multipoles together. 

7.1 Regular initial data 

If we assume that the invariants x, y are vanishing on the horizon then we get a nice hamiltonian 
system outside of the horizon. We can also assume that x and y are fixed and finite on the horizon. 
More precisely, we propose the following initial-boundary data: 

1. mon(x)=0, because of the singularity at r = 3m, moreover this charge indicates that we have 
chosen wrong parameter m in the background. 

2. dip(y) - weak internal angular momentum. 

3. The radiation data (x, X, y, Y) has to be finite on the horizon, moreover x, y should be controlled 
as a Dirichlet boundary condition. We assume also that = 0(1 /r) at spatial infinity. The center 

o 

of mass contained in dip[( A + 2) x] and linear momentum in dip(II) can be always "gauged out", 
see Appendix A.l, and those gauge transformations correspond to the infinitesimal translation and 
boost respectively. Performing those gauge transformations we pass to the "center mass rest frame" . 

4. The asymptotics of the invariants x, y at spatial infinity should guarantee finite hamiltonian and this 
can be achieved for standard asymptotics — 0(1 /r). Unfortunately, the standard asymptotics ~ 
does not guarantee finite angular momentum. More precisely, standard asymptotics at spatial infinity 



x,y = 0(l/r), x,y = 0(l/r ) gives logarithmically divergent integral in (7.4). The proposition 
by Christodolou and Kleinerman |Q, so-called S.A.F. condition, fits perfectly if we adopt it to the 
linearized theory. We propose the following asymptotics at spatial infinity 

x,y = 0(r- 3 / 2 ), x,y = 0(r- 5 / 2 ) 
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Remark Christodolou and Kleinerman assume that the full ADM data has better asymptotics (except 
conformal factor in the riemannian metric). Let us notice that if we assume that P kl = 0(r-^ 2 ) 
and hki = 0(r~ 3 ^ 2 ) then all boundary terms in the symplectic form (analyzed in Appendix B) are 
vanishing like - at spatial infinity. 

8 Conclusions 

We have shown a natural functionals which represent energy and angular momentum of the weak 
gravitational radiation on a Schwarzschild background. Particularly, we do not have to separate 
angular variables and the result is gauge-invariant. Moreover, the equations of motion for the gauge- 
invariant degrees of freedom correspond to the well-known Regge- Wheeler and Zerilli results. We give 
also a complete (together with the interpretation) analysis of the stationary solutions. 

After separation of the angular variables we should notice that the results presented in sections 4, 
6 and 7 are very close to the approach presented by Moncrief in ||. However, the invariants used by 
Moncrief in polar part are different and non-local. 



A Reduction of the initial data to the invariants 

It is convenient to introduce the following variables: 



II := 2rP 



3A 



\\A 



AP 3 



Q:=2h 3 + 2rh 3 A llA -rH^ 
The scalar constraint ( 4.17| ) takes the form 

^(rV^Q),3 + rV V* - J( A + 2)# - 
r 2 



° 6m , 
A+2 }tt = 



From vector constraint (4.15-4.16) we get 

6m 
r 



ry/v(y/vH) t 



A+2 - — ) rP 3A u + 2vr 2 S AB UB = 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



A.l Dipole polar part of the initial data 

In this subsection we consider only dipole parts of the variables and we denote by the same letter 
their dipole parts as the full o bject s i n the rest of the paper. 

1. From vector constraint (A.4), ( 4.15 ) and the definition ( |Al| ) we have 



rP 



3 A 



6m 



■(nMI), 



p\ = 



6m 2 



6 m 



■(V^n),, 



(A.5) 
(A.6) 
(A.7) 



This means that II contains the full information about dipole polar part of P 
2. From scalar constraint (A.3) we get 



h 3 3 



6m 
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Moreover, from the definition (A.2) we have 
2rh 3 A u A = Q + ^(r 2 ^Q),3 + rH. 3 



3 m 

This way we can see that Q and H contain the full information about dipole polar part of the linearized 
metric hki- 

3. The polar dipole gauge can be described by the following transformations: 



12mA 



-n 



12mA 



n + £° 



om om 



^ H + TV Q 
2 6m 



1 rv „ 

~H+—Q+e 

2 6m 



(A.8) 



and they show that we can always perform quasilocal gauge transformation in such a way that II, Q 
and H are vanishing. 

Moreover, from (4.8-4.12) we can check the evolution equations: 

o r 3 • 1 
6mA ^6^ 



^03 = 



12mA 
1 



(rH), 



12m 



Q 



6mA 



-n 



Finally we have shown that quasilocal gauge II = Q = H = gives vanishing dipole polar part of the 
full metric h„ 



A.2 Radiation polar part of the initial data in Regge- Wheeler gauge 

The special form of the metric proposed in [[lj will be called Regge- Wheeler gauge^. In polar part 
it gives the following gauge conditions: 

X AB \\ab = h A HA = h 3A W A = ; r\ AB UB — r\ AB UB + ( A + 2)£ A U 

r 2 h 0A U A —> r 2 h 0A W A + ACo + r 2 ^,| A ; r 2 h 3A ^ A — » r 2 h 3A W A + A 6 + r 2 ^ A U ), 3 (A.9) 

which allow to reconstruct polar part of the metric h^ v as follows: 
1. Equation 



2vK-'S A *\\ AB = x am \\ab - ( A + 2) W A 

follows directly from ( 4.1C ) and gives S AB \\ AB = 0. The variable II (defined by (AT)) contains the 
same information as X = 2r 2 S AB \\ AB + BJI. From ( |A.4|) we reconstruct P 3A \\ A , and (A.l) gives 
P 3 3 . Finally, we reconstruct S_ from vector constraint fl4.15|) . This way we have obtained four polar 
components of the ADM momentum P_ kl , namely S AB \\ AB , P_ 3A \\ A , P_ 3 3 and S_. The remaining two 
axial components are analyzed in subsection A. 4. 

!L Regge- Wheeler impose maximal number of gauge conditions: 1 axial gauge (A. 10) and 3 polar gauge conditions 

(A.9) in radiation part. They do not discuss mono-dipole part except some stationary solutions. Chandrasekhar in his 
book Q| does not impose axial gauge. Vishveshwara uses the same formalism as Regge- Wheeler plus axial gauge (A. 11) 
for the dipole and moreover in polar part he assumes dip(h[j + = 0. The monopole part is not discussed in the 
literature. 
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2. The spatial metric we reconstruct from the scalar constraint ( JA. 3| ) together with ( ]A.2j ) and the 

.. o 

observation that BQ = x + M A + 2)_ff. More precisely, we get a system of two equations for two 
missing components H_, h 3 ^ of the spatial metric (we need boundary data to solve them!). 

3. The lapse /ig we get from S AB \\ab which is given by (4.13). And finally the missing component 
h 03 of the shift we can calculate from the following equation 



2r 2 A~ 1 P 3A 



Ah 



r z h 



3 A \\A 



r 2 (h, 



OA 



which may be easily checked from (h9). The axial components of the spatial metric and axial part of 
the shift vector we analyze in subsection A. 4. 

The above analysis together with the results in subsection A. 4 show how to reconstruc t th e full 
i nitial data {P kl ,h kl ) together with lapse and shift h ok in the Regge- Wheeler gauge (A.9) and 
(A. 10) from the reduced initial data (x, y, X, Y). 



A. 3 Radiation polar part of the initial data in quasilocal gauge 

We do not like to impose any conditions directly on lapse and shift. It is more elegant to impose gauge 
conditions on the initial data only. The "time conservation laws" of the gauge conditions obtained 
from equations of motion give lapse and shift indirectly. For this purpose we propose the following 
quasilocal gauge conditions 

Q=K =n=o 

which allow to reconstruct radiation polar part of the four-metric in a quasilocal way. 



2 A 



X 1 O 

A n ^ A n -2 A 



6m \ 
r I 



e 

6 m 



A+2- — U° 



vQ 



A 



6m 

r 



H 



vQ 



A 



6m 
r 



H 



6m 

r 



Moreover, each component of the metric h depends quasilocally on the invariants. This can be 
shown as follows: 

1. From flA.2j ) and ( |A.3| ) we obtain metric components 
r 2 X AB \\AB =x, H = 



hi 



A 



6m 
r 



2. From flA.l| ), ( A.4 ) and ( 115 ) we get the ADM momentum 

6m 



9 2 cAB 
Lr °\\AB 



S = 



2vr A' 



A +2- 



A _1 X 



3. Moreover, 

rU - ArQ = 



2rP 3A u - 2« 



6m 
r 



X 



2 - — ) (rP 3 3 - 2Ah 03 ) 
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and 

uQ--(A+2 \H = A-^H - ( A + 2 — 



give the following components of the shift vector 



-l 



and finally the equation for II gives the lapse function 

n i ( ° 6m \ 1 / ° 
h° = A~ (A +2 ] ( Ax- 2vBx- 2rux 3 



The above quasiloca l for mulae for t he radiation polar part of the m etric h^ v allow to check the 
equations of motion ( [f.27j ) and ( 4.28 ) by inspection from the (ii^- 4.15 ). 

Why we prefer the quasilocal gauge? 
Because the initial data with compact support becomes more clear. The full (constrained) initial 
data with compact support gives the reduced initial data with compact support and the opposite is 
true only in quasilocal gauge. The data with compact support allows to avoid the boundary value 
problems. 

A. 4 Monopole part of the initial data 

The monopole part of the data seems to be not analyzed in the literature. We propose here the 
complete analysis of this simple "gap" . Let us prolongate the gauge condition H = which fixes the 
radial coordinate and for the monopole part let us assume 

mon(if) = 

This way mon(x) = 2ui3mon(/i 3 3) and the monopole part of the metric takes the form 

1 v° 
mon(H) = ; mon(/i 33 ) = -B -1 mon(x) = 2 — 

For time coordinate we propose the following gauge condition 
mon(P 3 3 ) = 

which is no longer quasilocal, it needs the boundary data at spatial infinity for the parabolic equation 
obtained from The similar situation (nonlocal reconstruction) we encounter during analysis of 

the dipole axial part (see the next subsection). 



From the vector constraint ( 4.15 ) and mon(P 3 3) = we obtain 
mon(S') = 

This way the whole ADM momentum is trivial in its monopole part. 
The trace of ( p~K| ): 



l^ = A -lp33 + ft ||A + ^ o3 



gives mon(/io3) = 0, and finally the lapse we get from monopole part of (4.11) 

r 2 ■ 

mon(2— P 3 3 ) = mon(-2ru/io, 3 + 2vhf) 

If we assume that h® vanishes at spatial infinity we obtain 
p o 

mon(/io) = — lnw 
m 

The infinity which we encounter in the monopole part of the invariant x at r = 3m suggests that 
we have to assume p° = and this way we exclude the possibility of manipulation with the mass 
parameter m in the background metric. 
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A. 5 Axial part of the initial data 

1. The ADM mo mentu m components P ZA ^ B e ab and S AC cb^a B we get from y and axial part of the 
vector constraint ( ff.l8| ). 

2. In radiation part (I > 1) we impose gauge condition 

X AC \\c B e A B = ; r 2 x Ac llCB s A B — r\ AC \ {CB e A B + ( A + 2)U\\ B e AB (A.10) 

which fixes £, A \\ B s AB quasilocally. The component h 3 A\\B £AB is obviously contained in Y. Moreover, 
from ( 4.1 0| ) we have 



r X ||cs£a — ||cs£a + ( A + Z)Ii 0A \\b£ 



which gives h A\\B£ AB = 0. This way we have shown how to reconstruct axial part of the metric h 
from the invariants fy, Y). 

3. The dipole part of the metric (1 = 1) can be fixed by the gauge condition 

dip(h 3A \\Be AB ) = ; dip(^3A||B£ AB ) — » dip(h 3 A\\ B e AB ) + [dip(U\\ B z AB )] , 3 (A.ll) 

which gives "parabolic" equation for the angular gauge transformation dip(^A\\B£ AB ) with the bound- 
ary value at spatial infinity. Moreover, stationary solution in dip(y) appears in h A\\B£ AB ■ More 
precisely, from "time conservation law of gauge condition" 

= r 2 dip(h 3A \\B£ AB ) = r 2 dip(h 0A \\B£ AB ),3 + dip(y) 

4s 



we obtain dipole axial part of the shift dip(/io J 4||s£ ) = — , and in particular for s = scos9 we have 

,2 , 



2s sin 9 

Iiq^ = which has been proposed by Vishveshwara (eq. 5.2 in ||). 

B Reduction of the symplectic form 

Let (p kl ,qki) denotes the Cauchy data on a hypersurface E. The (2+l)-splitting of the tensor qki 

2 

gives the following components on a sphere: Q:= r\ AB qAB, <?33 ~ scalars on S 2 , q 3 A - vector and 

2 

1 ab ■= qAB — 2VAB 1 - symmetric traceless tensor on S . Similarly, we can decompose the tensor 
density p kl . On each sphere S(r) we can manipulate as follows 

/ p kl qki = [ p 33 q33 + 2 P 3A q 3 A + lm + P AB ^AB = 

Jv Jv * 

= f^q 33 - 2(rp* A u ) A~\rq 3A llA ) 2(rp^ B e AB ) k~' \rq 3A \\ B e AB ) +\% + 

+2 f (r 2 e AC Pa B \\bc) ( A + 2)-\r 2 e AC q A B \\Bc) + 
Jv 

+2 / (r 2 P AB \\ AB ) A _1 ( A + 2)-\r 2 °q AB \\ AB ) 
Jv 

where we have used the following identities on a sphere 

- / n A v A = (rTT A \\ A ) A~\rv A llA ) + (rn A tt B e AB ) A ^ (rv A \\ B e AB ) (B.l) 

Js(r) 

and similarly for the traceless tensors we have 

/ 7r AB v ab = 2 f (r 2 e AC n A B l{BC ) A ~\ A + 2)-\r 2 s AC °v A B \\bc) + 

JS(r) Js(r) 
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+ 2 (r 2 n AB llAB ) A ( A + 2)-\r 2 v AB ]]AB ) 

JS(r) 

The axial part of the quadratic form J v p kl qki we define as follows: 

axial part = -2 / (rp 3AHB e AB ) A (rq m \ B s AB )+ 
Jv 



(B.2) 



+2 / (r 2 e AC Va B \\bc) A ( A + 2Y\r\ AC °q A B 



\\BC) 



The remainder we define as a polar part. Moreover, from the axial part of the vector constraint (4.18): 



[r p " £AB),3+re Pa \\bc — u 
we obtain 

axial part = -2 / (rp 3A ^ B e AB ) A (rq 3 A||s£ AB )+ 



-2 / (rV A|l ^AB), 3 A ( A + 2)~ {r e 9a ||sc) ; 
= -2 [ {rV A ^ B SAB) {k+2)-\r*e AC q A B \\Bc) 

JdV 



-2 / (rV A|l ^AB)A feA|| S ^-(A+2)- i (r^g^|| BC ),3] 

We can see the invariants in the volume term if we write dipole part separately 

dipole axial part = -2 / dip(r 2 p 3AllB e A b) k^ 1 {q iA \\B£ AB ) 
Jv 

and finally the radiation axial part contains gauge-invariants in the volume term 

radiation axial part = -2 f {r 2 p 3A ^ B e A B) A ( A + 2y 1 (r 2 e AC q A B \\bc) + 

JdV 



-2 / (r 2 p 3AllB SAB)A (A+2)- 1 [(A+2)q 3A]]B e AB -(r 2 e A ^ q A B \\BcU\ 
Jv 

For the polar part we can use the rest of the vector constraints 

-^(Vvp 3 3 ),3 + rp 3 A llA - P= 
y/v 

(r 2 p 3A HA),3 + (r 2 P AB U B) + ^kp=0 
and we can reduce partially polar part as follows 

polar part = j^m* ~ 2(rp 3A u ) A + \ (^=(V^p 3 3),3 + rp 3A \\ A ) 5 



(r 2 p 3A \\ A ),3 + \k (^=(V^p 3 3),3 + rp 3 A U 



(B.3) 



(B.4) 



(B.5) 

(B.6) 
(B.7) 



A (A+2)-t(r 2 q AB llAB 



rp 3 3 



av 



1 2 o 



9-(A+2)- 1 (r 2 <Z AB || AB 



-2 I (rW)A (A+2)-\r 2 °q AB llAB ) + 



av 
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rp3A' 



±q+2rvA \A+2)- 1 (r 2 q AB llAB ), 3 -(A+2)- 1 (r 2 °q AB llAB )-2A \rq 3A ]]A ) 



pS 



Q . o 1 / 7 ° 

q 3 3 + V^{A+2)- 1 —q 



'AB 



\\AB I > 3 



v I r 2 



9 ,3 



(B. 



The above calculation shows that we should consider mono-dipole part separately. The monopole part 
is very simple 



mon( / p kl qu) = [ \p 33 B 1 mon(C) + - / rp 3 3 mon(q) 

JV JV Z 1 JdV 

where invariant C is defined as follows 



(B.9) 



(:=B 



2q 33 + 2rq 3A n A -rvq, 3 



r 2 q AB UB -\{A+2) 2 q 



Using (A. 6), (A. 7), (A. 5) and integrating by parts we can also rewrite (from the beginning) dipole 
part in the following way 

r — 1 1 22 

dipole polar part = / p 33 q 33 - 2(rp 3A \\ A ) A {rq 3A \\ A ) + - pq= 



V 



^n,, 3 -in 



933 - 2 "7Z~"(v / ^n) ;3 A (rq 3A \\ A ) + - 



J dv 12m 6m \ 



6m 



iA \\A + q 33 - \vr g, 3 



6 m 



2 

q= 











6m 


r 2 ^( 



1 2 



}=/ ^(V^+^iiq 

,J J 9V 12m 12m 



rvU A -1 A + 2 

av \ r 



6m \ 1 



o , 2 



(B.10) 



where here II d enotes only its dipole part and II itself is defined by (A.l). We h ave also used scalar 
constraint ( |A.3| ), and finally the dipole polar part takes its boundary form ( B.10|). 

The radiation pol ar pa rt can be also reduced if we use scalar constraint ( 4.17 ) in two equivalent 
forms (motivated by (B.8)) 



v / r 2 



«s° - -V -5= 9 , 3 



yfi^—q AB l[AB ) ,3 + (A +2) 
~ A( A + 2) q +2rv (r 2 °q AB \\ AB ) ,3 - Ar 2 q AB \ ]AB - 2( A + 2) (rq 3A u ) 



(B.ll) 



= 2v^ f — Cj , 3 + 2«£C - ( A + 2)B- 1 ( 
Inserting ( B.ll ), ( B.12| ) into (B^) and integrating by parts we obtain 



radiation polar part = / p 3 3 



rp 3 A' 



(B.12) 



<Z3 3 + Vv( A + 2)- 1 — q AB \\ AB ,3 - -V "7= « » 3 



v r 2 



2 V> 



i I +2r, A~ 1 ( A + 2)- 1 (r 2 g^,,^) ,3 - ( A + 2)- 1 (r 2 



-2 A (rg 3 V) 
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+ / rp° 3 

JdV 



1 2 o / 

-9-(A+2)- 1 ( 



(rW) A _1 (A+2)- 1 (r 2 ^ B ||AB ) 



/ rp 3 3 i 3 -rU A ( A + 2)- 1 (r 2 9 AB ||ab) + / p\( A + 2)" 
Jav ^ v y 



rp 3A ,uA (A+2)- 1 



ar 



p 3 + 2 A (rp 3 



2y^ ( ^=C ) , 3 + 2BC - ( A + 2) U - 1 £T 1 C 



(y/vR) 3 +Sn-( A +2)6- 1 w- 1 rp 3A || j4 



1 2 



av ^ 



( A + 2)-\ + rp%- q -rU A ( A + 2)" 1 (r 2 9 AS ||ab) + 
A (A+2)- x C = 

/ [2r 2 P AB wab+BU) A ( A+2)- 1 C(B.13) 



9+rnA ( A+2)" 1 (c-r 2 



where we used equality 



2r 2 + -± 3 + ( A + 2)ZT 1 «- V A 



1 1-4 







which is a simple consequence of the vector constraint (A. 4), ( JB.6| ) and (B.7). It is easy to check that 
the volume terms in the final form of (B.5) and (B.13) contain demanded gauge-invariant result. 
Let us summarize the result we have proved 



P Qkl = monopolc part + dipolc part + radiation part 

monopole part is given by ( |B.9| ), the dipole part is a sum of (B.4) and ( B.1C| ): 

dipolc part = dipole axial part + dipole polar part 
and radiation part is a sum of ( |B.5| ) and ( B.13| ). 

" A _1 ( A +2)" X C+ 



radiation part in V = 



2r 2 P AB \\ AB + BU 



-2 / (rV AllB e AB ) A" 1 q 3A \\BS AB -(A+2)-i(rh Ac qA B llBC ), 3 
Jv L 

f ° i ( 6m \ 1 o , A 2 
boundary terms = J rvll A ( A+2 \ Q-r 2 p ||A A 9 + 



-2 / (rV A|l ^ S ) A ( A + 2)- > 2 e^ 
./ay 



C (2+l)-decomposition of the gauge for lapse and shift 



The gauge transformation ( |2.18| ) splits according to the (2+l)-decomposition and we obtain the fol- 
lowing gauge transformation law which acts on lapse and shift: 



2ttl 

hoo — > h 00 + 2£ — ^"£3 

hoA — > hoA + Co, A + £a,o 
ho3 -> h 03 +i 3 - v£° 3 



(C.l) 
(C.2) 
(C.3) 
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